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Abstract 

Minimal surfaces with isothermal parameters admitting Bezier representation were 
studied by Cosm and Monterde. They showed that, up to an afhne transformation, 
the Enneper surface is the only bi-cubic isothermal minimal surface. Here we study 
bi-quartic isothermal minimal surfaces and establish the general form of their gen¬ 
erating functions in the Weierstrass representation formula. We apply an approach 
proposed by Ganchev to compute the normal curvature and show that, in contrast 
to the bi-cubic case, there is a variety of bi-quartic isothermal minimal surfaces. 

Based on the Bezier representation we establish some geometric properties of the 
bi-quartic harmonic surfaces. Numerical experiments are visualized and presented 
to illustrate and support our results. 

Key words: minimal surface, isothermal parameters, Bezier surface 

1 Introduction 

Minimal surfaces have recently become subject of intensive study in physical and bi¬ 
ological sciences, e.g. materials science and molecular engineering which is due to its 
area minimizing property. They are used in modeling physical phenomena as soap dims, 
block copolymers, protein folding, solar cells, nanoporous membranes, etc. Minimal sur¬ 
faces found applications also in architecture, CAGD, and computer graphics where Bezier 
polynomials and splines are widely used to efficiently describing, representing and visu¬ 
alizing 3D objects. Hence it is important to know minimal surfaces in polynomial form 
of lower degrees. Bi-cubic polynomial minimal surfaces are studied in |1]. Polynomial 
surfaces of degree 5 and 6 are studied in [6] and where some interesting surfaces are 
described and their properties are examined. Examples of polynomial minimal surfaces 
of arbitrary degree are presented in [8]. 

In section 3 we specify the result of Cosm and Monterde [1] concerning bi-cubic poly¬ 
nomial minimal surfaces. We note that their proof that these surfaces coincide up to 
affine transformation with the classical Enneper surface concerns the case of surfaces in 
isothermal parameters. In section 4 we consider an analogous problem for polynomial 
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surfaces defined by charts x(m, v) of degree 4 on both u and v. It turns out that they are 
more various than the bi-cubic ones, so they may be more useful in computer graphics. In 
this paper we use a new approach to minimal surfaces proposed by Ganchev [2] as well as 
the method from |3] to obtain a parametrization of the surface in canonical parameters. 

In section 5 we consider bi-quartic harmonic Bezier surfaces. We show that for a 
special choice of nine boundary control points the corresponding harmonic Bezier surface 
is uniquely determined and is symmetric with respect to one of the coordinate planes 
Oxy, Oxz, and Oyz. Based on the Bezier representation we apply computer modeling 
and visualization tools to illustrate and support our results. 

2 Preliminaries 

Let S' be a regular surface. Then S is locally dehned by a chart 

X = x(m, v) {u, v) e U C . 

As usual we denote by x^, Xy, x„„,... the partial derivatives of the vector function x(u, v). 
Then the coefficients of the hrst fundamental form are given by the inner products 

E = xl , F = XyXy , G = xl 


and the unit normal is 

^ _ x„ X x„ 

|Xn X X^l 

Then the coefficients of the second fundamental form are dehned by 

E U Xyy^ AL U Xyy^ TV U Xyy . 

The Gauss curvature K and the mean curvature Lf of S' are given by 

_LN -M^ ^ _EN - 2FM + GL 

~ EG-F^' ~ 2{EG - F2) ’ 

respectively. Note that the Gauss curvature and the mean curvature of a surface do not 
depend on the chart. The surface S is said to be minimal if its mean curvature vanishes 
identically. In this case the Gauss curvature is negative and the normal curvature of S is 
the function u = see [2]. 

We say that the chart x(m, v) is isothermal or that the parameters (u, v) are isothermal 
if ST = G, F = 0. It is always possible to change the parameters (u, v) so that the resulting 
chart be isothermal. We note however that this change of the parameters to isothermal 
ones is in general nonlinear. 
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When the chart is isothermal it is possible to use complex functions to investigate it. 
We shall explain briefly this. Namely let f{z) and g{z) be two holomorphic functions 
(actually sometimes they are taken meromorphic). Dehne the Weierstrass complex curve 
by 

^( 2 :) = ^ 9 \z)),y{z){l +g‘^{z))J{z)g{z)^ dz . (1) 

Then T(; 2 ) is a minimal curve, i.e. (\h'(z))^ = 0, and its real and imaginary parts 
x(m, n) = i?e ^( 2 ;) and y(M, n) = Jm T(z) 


are minimal charts. Moreover, they are isothermal and are harmonic functions (i.e. Ax = 
0, Ay = 0, where A is the Laplace operator) as the real and complex part of a holomorphic 
function. Conversely, every minimal surface can be dehned at least locally in this way. 
Of course a minimal surface can be generated by the Weierstrass formula with different 
pairs of complex functions f{z), g{z). 

It is easy to see that the coefficients of the hrst fundamental form of a chart dehned 
via the Weierstrass formula with functions f{z), g{z) are given by 

B = G = i|/pl + |jPy F = 0. (2) 

The normal curvature is computed to be 

^ Ma'l .ON 

m + \9?r 


see [3], Theorem 22.33. 

Recently Ganchev [2] has proposed a new approach to minimal surfaces. Briehy speak¬ 
ing he introduces special parameters called canonical principal parameters. A chart with 
such parameters is isothermal. Moreover, the coefficients of the two fundamental forms 
are given by 

F = 0, G=i, 

L = l, M = 0, N = -1. 

His idea leads to the fact that the real part of the minimal curve 


<h(tc) 


rni-g^{w) il + g^{w) 9{w)\ , 
Ao V2 ~9» ’2 ~g'{w) '~g'{w)J ^ 


( 4 ) 


is a minimal surface in canonical principal parameters. Note that this is the Weierstrass 
formula with f{z) = —l/g'{z), g{z) = g{z). 
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We shall use also the following theorems: 

Theorem A. [2] If a surface is parametrized with canonical principal parameters, 
then the normal curvature satishes the equation 


A In z/ + 2z/ = 0. 


(5) 


Conversely, for any solution u{u, v) of equation ([^ there exists a unique (up to position 
in the space) minimal surface with normal curvature u{u,v), where {u,v) are canonical 
principal parameters. □ 

Theorem B. |1] Let the minimal surface S be dehned by the real part of ([^. Any 
solution of the differential equation 


{Aw)r 


1 

f{z{w))g'{z{w)) 


( 6 ) 


dehnes a change of the isothermal parameters of S to canonical principal parameters. 
Moreover, the function g{z) that dehnes S via the Ganchev’s formula (|^ is given by 
g{w) = g{z{w)). □ 

The canonical principal parameters [u, v) are determined uniquely up to the changes 


u eu + a, ^ ^ 

V = ev + 0 ^ 


3 Bi-cubic minimal surfaces 


By investigating minimal Bezier surfaces Cosin and Monterde [T] formulate that any bi¬ 
cubic minimal surface dehned by 

/ 3 3 3 \ 

x(u,v) = j bijU^v^, CijU^v^ I 

\i,j=0 i,j=0 i,j=0 / 

is, up to affine transformation in the space, actually an affine reparametrization of the 
classical Enneper surface 

enneper(M, t') = - I m —v -1- —— — v‘^ \ . 

Z o o J 

(This chart of the Enneper surface is obtained from the Weierstrass formula with f{z) = 1, 
g{z) = z.) We note that their proof actually refers to bi-cubic isothermal minimal charts. 
Indeed, as we mentioned in Section the change of parameters to isothermal ones is in 
general nonlinear. Below we give a simple example of a bi-cubic minimal chart that can 
not be transformed by an affine transformation into an isothermal one. 
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Example. Consider the bi-cubic chart 


If o 2 2 2 2 

x(m, v) = - \uv ---h uv , —V + —— U V ,U V — V 

Z o o 

It can be shown by direct computation that this chart dehnes a minimal surface. Of 
course we can simply remark that this is a reparametrization of enneper(M,n) with u 
replaced by uv, so the mean curvature vanishes identically. 

Let us make an affine transformation of the parameters {u,v)-. 


u = aiu + hiv + Cl 

V = a2U + b2V + C2 


with nonzero Jacobian, i.e. 

J = aib2 — 0261 7 ^ 0. ( 8 ) 

We shall try to determine the coefficients ai,bi,Ci, so that the chart 

x(-u, v) = x(ai'U + biv + Cl, a 2 U + 62 ^ + C 2 ) 

be isothermal. Actually we shall see what follows only from F = 0 . A direct computation 
shows that F = = \Fi ■ F 2 , where 

Fi= + (c2 + a2U + b2vY{l + (ci + aiu + 6ih)^)j , 

F2 = albiu{ 2 aiu + biv + ci) + 01(62^ + C2)(6iC2 + b 2 {aiu + 2 biv + ci)) 

+02 ^6iC2(3ai'U + biv + Ci) + 62(1 + c? + SoiCih + 2a\u^ + 36i(ci + 2 ai'u)h + 2b\v‘^^ 

Since Fi is positive, the vanishing of F implies F 2 = 0. Hence the coefficients in F 2 must 
be zero. In particular the coefficients of , uv are 

aia2(a2&i + 0162) = 0, 
bib2{a2bi + 0162 ) = 0 , 

( 02^1 + 01 ^ 2 )^ A 40102^1^2 = 0 - 

These equations imply immediately 

fllQ' 2 &l &2 = 0, 02^1 + 0 - 1^2 = 0. (9) 

Let e.g. oi = 0. From (|^ it follows 02^1 = 0, which contradicts ([^. So it is impossible 
to make an affine transformation of the parameters in Q to obtain an isothermal chart. 

Remark. In view of the above notes the problem of existing bi-cubic minimal surfaces 
different from the Enneper one is still open. More generally it will be interesting to obtain 
a method for finding polynomial minimal non-isothermal charts. 
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4 Bi-quartic minimal surfaces in isothermal parame¬ 
ters 


In this section we examine minimal snrfaces represented by isothermal polynomial charts 
of degree 4 in both n, v. We may expect that there exists more than one snch snrface, 
but it is interesting to know “how many” are there. 

So consider the chart 

4 

x(M,n) = ^ VijuV, 

i,j=0 

where Vjj = {aij, bij, Cij) are vectors in Using F = 0 and looking on its coefficient of 
we obtain V 44 = 0 . Analogously we derive consecutively V 43 = 0 , V 34 = 0 , V 42 = 0 , 

V24 = 0, V41 = 0, vi4 = 0, V33 = 0, V32 = 0, V23 = 0. 

It is known that any minimal isothermal chart is harmonic, see e.g. [3]. In our case 
this implies 

V02 = ~V20, V21 = — 3 Vo 3 , V22 = —6V40, 

V12 = —3 v 30 , Vi3 = —V31, V04 = V40. 

Substituting these in F and looking on the coefficients of u® and u^v we obtain also 
V 31 V 40 = 0, — 16 v|q = 0. If V 40 = o the chart is not of degree 4. So we assume 

V 40 7 ^ o. Up to position in space and symmetry we may take 


V4o = (p,0,0), V31 = (0,4p,0), where p 7 ^ 0 . 


Now the coefficients of u® and u^v in F give 693 + O 30 = 0 and 093 — 639 = 0. Using 
this we can calculate the derivatives and x„. Let the functions f{z) and g{z) give 
the Weierstrass representation of the surface. Denote by (0i,02,03) the derivative of 4/. 
Then T' = (0i, 02, ^ 3 ) = — ix^. In our case a direct computation shows 

01 = -m9i + ai9 + (m + iv) (^-iau + 2a2o + {u + iv) (3039 + 3*639 + 4 p(m + in)) j, 

02 = -ihoi + 619 - i{u + iv) 1^611 + 2*620 + (n + iv)[ 3 a 3 o + 3*639 + 4p(** + *n)) 

03 = —icoi + Ci9 + (—*Cii + 2 c 29 + 3(*C93 + C39)(** + *n)) (** + iv). 

On the other hand the Weierstrass formula implies easily 

/(^)=01-*02, ^(^) = 

01 - *02 

Hence we derive 

f{z) = —*091 — 691 + 019 ~ *619 + (— *Oii — 611 + 2o29 ~ 2 * 629 ) 2 ;, 

, _ C91 + *Ci9 + (Cii + 2 *C 29 + 3 ( —C93 + ic 3 o)z)z 

Ool — *691 + *Oi9 + 619 + (Oii — *611 + 2*029 T 2629)^ 
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Consequently we have obtained that for some complex constants A and B 

f(z) = Az + B, g(z) = 

where P2{,z) is a polynomial of degree at most 2. Suppose A = 0 , i.e. f{z) is a constant. 
Then the derivative 

{K 02 , 03 ) = (^))> 

is of degree 2 or 4 , so the chart x(m, v) is of degree 3 or 5 , which is not our case. So 
A ^ 0 . Since 0 i = ^f{z){l — g^{z)) is a polynomial then Az + B divides P2{z). Hence 
g{z) = Cz + D, where C 7^ 0 . We have proved the following 

Theorem 1 Any bi-quartic parametric polynomial minimal surface in isothermal param¬ 
eters is generated by the Weierstrass formula with the functions 


f{z) = Az + H, g{z) = Cz + D, where A 0 , C 7^ 0 . 


Further, we are interested which of the functions in Theorem generate different surfaces. 
Denote by xo(m, n) the chart defined as the real part of the Weierstrass minimal curve 
with functions f{z) = z, g{z) = z and the corresponding surface by So- Denote also by 
x(M,n) the chart with generating functions f{z) = Az, g{z) = Cz for arbitrary nonzero 
complex numbers A, C, and the corresponding surface by S. Using ([^ and ([^ we can 
see that the nonzero coefficients of the hrst fundamental form and the normal curvature 
of xo(M,n) are respectively 

4 

+ v‘^{l + 

Obviously xq(u,v) is not in canonical principal parameters. We want to change the 
parameters (u, v) to canonical principal ones. Equation ([^ has a solution 

z = (3/2)2/3 (^w)2/3_ 


Eq = Go = -{m + n^)(l + , and hence z/q = 




So according to Theorem B we change the variable 2; by ( 3 / 2 )^/^(z zY^"^■ Now the functions 


Rz) = t 



g{z) = 



(iz) 


2/3 


generate a chart xo(n,n) in canonical principal parameters and 


Vo = 


(^2 + ^2)1/3 (l + (|)U3 ^^2 + ^2y/3y 
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Analogously x(m, v) is not in canonical principal parameters. Using again Theorem B 
and changing the complex variable 2; by 


3 \ 


2/3 


( 


i z 


2/3 


) 


we obtain a corresponding chart x(m, v) in canonical principal parameters. According to 
(|^ its normal curvature is 



The last formula implies that this is also the normal curvature (in canonical principal 
parameters) of the surface Si generating via the Weierstrass formula by the functions 


hi.^) 



gi{z) = z. 


According to Theorem A the surfaces S and Si coincide up to position in the space. On the 
other hand, the Weierstrass formula implies that the surfaces Sq and Si are homothetic. 
So for any nonzero complex numbers A, B the surface S is, up to position in the space, 
homothetic to Sq. Surfaces of type S for different values of A and C are shown in Figure 



Generating functions f(z)=10z, g(z)=z Generating functions f(z)=z, g(z)=10z 

Figure 1 : Surfaces of type S for different values of A and C 

Now we are interested whether the functions f{z) = z, g{z) = z + a + 16, a, 6 G M, 
dehne a surface which is really different from Sq. The chart :si{u,v) generated by these 








functions is not in canonical principal parameters so according to Theorem B we change 
the complex variable by ( 3 / 2 )^/^(i Then the functions 


f{z) = i 



1/3 

g{z) 



2/3 

{i + a + ib 


define a chart x(m, v) in canonical principal parameters. Its normal curvature is 


u = 


4 ( 1 ) 


2 \ 2/3 


\/u‘^ + '^^(1 + B{u, v)B{u, v)y' 


where B{u,v) is 


/ 3 \ 

B{u,v) = a + ib + l-j {iu — vY^^. 


Comparing these functions for different values of (a, b) we can say that the resulting 
surfaces are different. In Figure]^ are shown parts of the surface S'o, obtained for (a, b) = 
( 0 , 0 ) (left), the surface obtained for (a, 6) = ( 0 . 5 , 0 ) (center), and the surface obtained 
for (a, 6) = (1,0) (right). 



f{z) = z, g{z) = z f{z) = z, g{z) = z+ 1/2 f{z)=z,g{z) = z + l 

Figure 2 : Comparison of bi-quartic minimal surfaces for generating functions f{z) = z 
and g{z) = z + a + ib 


5 Bi-quartic harmonic Bezier surfaces 

We consider bi-quartic tensor product Bezier surface defined by 

4 4 

x(m, -y) = h,jBf{u)Bl{v), (10) 

i=0 j=0 
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where bj^, i, j = 0,... ,4 are the control points of x(m, n), and Bf{u) are the Bernstein 
polynomials of degree 4 dehned for 0 < m < 1 by 


Bt{n) :=( 


4—i 



lo, 


for i = 0,..., 4, 
otherwise. 


Recall that if x(m, v) is in isothermal parameters then x(m, v) is a minimal snrface if and 
only if x(m, v) is a harmonic surface, i.e. Ax = 0. For a harmonic Bezier surface Monterde 
[S] has proved that if we know the control points on two opposite boundaries except one 
corner point, e.g. nine points {boj}|=o {bi 4 }f= 0 ) then the remaining sixteen control 
points are fully determined. The prooQis based on the harmonic condition Ax = 0 which 
leads to a linear system that has a unique solution. Here we assume that nine control 


b40 O O O b44 

bso o O O b34 

b 20 o o o o 

bio o o o bi 4 

boo o o o bo4 

Figure 3: Input control points that fully determined bi-quartic harmonic Bezier surface 


points {bio}i=o! {bi 4 }i=o;i ^2 shown on Figure are given. Note that they differ from 
those used in [5] . In Lemma below we give expressions for the remaining control points 
through the given points. Then in Proposition!^ we prove that in the case where the given 
points are symmetric with respect to any of the coordinate planes, then the corresponding 
harmonic Bezier surface is symmetric with respect to the same coordinate plane. 


Lemma 1 Let nine control points b^o, i = 0, ...,4; bj 4 , i = 0,1,3,4 be given. A bi- 
quartic Bezier surface (10) is harmonic if and only if the remaining sixteen control points 
satisfy 


^Monterde’s proof is made for a surface of degree (n,n), where n G N is even. Here we consider the 
case n = 4. 
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boi — (25boo + 8bo4 — 40bio — 8bi4 + 36b2o — ISbao + 4b34 + 4b4o — b44)/12 

b4i = (4boo ~ bo4 — 16bio + 4bi4 + 36b2o ~ 40b3o — 8b34 + 25b4o + 8b44)/12 

bii = (ITboo + 7bo4 — 14bio — 4bi4 + 18b2o — 8b3o + 2b34 + 5b4o + b44)/24 

bsi = (5boo + bo4 — 8 bio + 2bi4 + 18b2o — 14b3o — 4b34 + 17b4o + 7b44)/24 

bo 2 = (13boo + 8bo4 — 28bio — 8bi4 + 30b2o ~ 16b3o + 4b34 + 4b4o — b44)/6 

b42 = (4boo ~ bo4 — 16bio + 4bi4 + 30b2o ~ 28b3o — 8b34 + 13b4o + 8b44)/6 

bi 2 = (llboo + 7bo4 — 20 bio — 4bi4 + 24b2o — 14b3o + 2b34 + 5b4o + b44)/12 

b32 = (5boo + bo4 — 14bio + 2bi4 + 24b2o — 20b3o — 4b34 + llb4o + 7b44)/12 

bo3 = (14boo + 19bo4 — 32bio — 16bi4 + 36b2o — 20b3o + 8b34 + 5b4Q — 2b44)/12 

b43 = (5boo ~ 2bo4 — 20bio + 8bi4 + 36b2o ~ 32b3o — 16b34 + 14b4o + 19b44)/12 

bi3 = (5boo + 7bo4 — 8 bio — bi4 + 9b2o — 5b3o + 2b34 + 2b4o + b44)/12 

bss = (2boo + bo4 — 5bio — 2bi4 + 9b2o — 8b3o — b34 + 5b4o + 7b44)/12 

b2i = (3boo + bo4 — 4bio + 8b2o — 4b3o + 3b4o + b44)/8 

b22 = (7boo + 3bo4 — 16bio + 24b2o ~ 16b3o + 7b4o + 3b44)/12 

b23 = (7boo + 5bo4 — 16bio + 4bi4 + 24b2o — 16b3o + 4b34 + 7b4o + 5b44)/24 

b24 = (boo ~ bo 4 — 4bio + 4 bi 4 + 6 b 2 o — 4 b 3 o + 4 b 34 + b 4 o — b 44 )/ 6 . 

( 11 ) 

Proof. It follows straightforward using the corresponding linear system from [S]. □ 

Proposition 1 Let the given points bjo, i = 0,...,4; bj 4 , i = 0,1,3,4 be symmetric 
with respect to some of the coordinate planes Oxy, Oxz, or Oyz. Then the corresponding 
harmonic Bezier surface defined by Lemma\^is symmetric with respect to the same plane. 

Proof. Let hij = hij{xij, yij, Zij) and assume that the given control points are symmetric 
with respect to the plane Oxy, i.e. bofc and bi^ are symmetric points to h^k and bo^, 
respectively, /c = 0,4, and b 2 o lies on Oxy. Then we have 


Xok — a^4fc, Xik — Xsk, yok — UAk, Ulk — Usk, ^Ok — —^Ak^ ^Ik — —^^k 


( 12 ) 


for k = 0,4, and ^20 = 0. To show that the harmonic Bezier surface defined by Lemma 
is symmetric with respect to Oxy it suffices to establish that its control points are 
symmetric with respect to Oxy. We need to establish that ( 12 ) holds for k = 1,2,3 and 


Z 2 j = 0 ioT j = 1,..., 4. Next we verify that Xqi = Xn, y^i = 1 / 41 , Zqi = — 2 : 41 , and 2:21 = 0. 
The analogous relations for the remaining control points follow in a similar way. 

From 0 and ([T^ we have 


xoi = (25x00 + 8 x 04 - 40xio - 8 x 14 + 36x20 - 16 x 30 + 4x34 + 4x4o - 2 ) 44)/12 

= (29xoo + 7 xo 4 - 56xio - 4 xi 4 + 36x2o)/12 = X 41 . 

Analogous relation holds for yQi and 1 / 41 . For the third coordinates 2:01 and 2:41 we obtain 


Zqi — (21x00 + 9xo 4 — 24xio ~ 12xi4)/12 — —X 41 . 
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,'512 176 128 \ 

1 3 ’ 3 ’ 3 1 

fl28 424 152 '1 

1 3 ’ 3 ’ 3 t 

(-128, f,f) 

f 128 296 104 '1 

V 3 ’ 3 ’ 3 / 

,'256 80 128 

1 3 ’ 3 ’ 3 

64 536 88\ 

3 ’ 3 ’ 3 

f 224 44 32 \ 

1 3 ’ 3 ’ 3 

(0,0,1) 

/160 20 32 \ 

V 3 ’ 3 ’ 3 

,'64 280 401 

V 3 ’ 3 ’ 3 

/ 512 r| 160 \ 

(0,0,-I) 

(-lf,0,0) 

(0,0,-I) 

(-f^,0,32) 

f 64 536 88 \ 

1 3 ’ 3 ’ 3 

f 224 44 32 

1 3 ’ 3 ’ 3 

(0,0,1) 

f 160 20 32 '1 

1 3 ’ 3 ’ 3 t 

(64 280 40\ 

V 3 ’ 3 ’ 3 

,'512 176 1281 

V 3 ’ 3 ’ 3 t 

,'128 424 152 \ 

V 3 ’ 3 ’ 3 

(-128,-If, f) 

( 128 296 104 \ 

V 3 ’ 3 ’ 31 

,'256 80 1281 
V 3 ’ 3 ’ 31 


Table 1: Control points of a harmonic bi-quartic Bezier surface that are symmetric with 
respect to Oxz 

It remains to show that 2:21 = 0. We have 

Z2I = ( 3^00 + ^04 — 42:10 + 8^20 ~ 42:30 + 3^40 + ^44)/8 

= (3(^00 + Z40) + (^04 + ^44) — 4(2:10 + ^ 3 o) + 82:20)78 = 0 . 

The case where the nine given points are symmetric with respect to the other coordinate 
planes is treated analogously. □ 

A bi-quartic harmonic Bezier surface which is symmetric with respect to Oxz is shown 
from two different viewpoints in Figure]^ Its control points are presented in TableWe 
note that they are obtained from the minimal bi-quartic Bezier surface with generating 
functions f{z) = z, g{z) = z — 1. Hence, the surface in Figure is harmonic minimal 
Bezier surface. 


6 Conclusions and Future Work 

In this paper we characterize all bi-quartic parametric polynomial minimal surfaces by 
their generating functions using the Weierstrass formula. We also consider the bi-quartic 
harmonic Bezier surfaces and establish their symmetry with respect to any of the coordi¬ 
nate planes. We present numerical experiments and give examples. A possible direction 
for future work is to extend our results for minimal surfaces of higher degrees. 
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